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Abstract 

In this paper limiting distribution functions of field and density fluctuations are ex- 
plicitly and rigorously computed for the different phases of the Bose gas. Several 
Gaussian and non-Gaussian distribution functions are obtained and the dependence 
on boundary conditions is explicitly derived. The model under consideration is the 
free Bose gas subjected to attractive boundary conditions, such boundary conditions 
yield a gap in the spectrum. The presence of a spectral gap and the method of the 
coupled thermodynamic limits are the new aspects of this work, leading to new scal- 
ing exponents and new fluctuation distribution functions. 
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1 Introduction 



Normal and critical fluctuations in the ideal Bose gas with Dirichlet or periodic boundary 
conditions are explicitly studied at different levels of rigour in |I|-[TO[|. Condensation occurs 



only at dimensions u > 3. It turns out that the behaviour of the energy gap of the finite 
volume spectrum as a function of the volume is determining for the degree of abnormality 
of the fluctuations in the condensate regime. Typical is that the limit spectrum has no 



energy gap. The ideal Bose gas with attractive boundary conditions 0,|T2[ on the other 
hand is quite different in nature. Condensation is possible in all space dimensions, and the 
spectrum has a finite energy gap in the thermodynamic limit. In this paper we analyse 
the nature of field and density fluctuations for the Bose gas with attractive boundary 
conditions. Field fluctuations are centred observables linear in the Bose-field operators 
and density-fluctuations are quadratic in the Bose-fields. We derive rigorously the exact 
form of the limiting characteristic functions, i.e. we study the limits: 

lim ude'^^'') , for t G R, (1.1) 

where L indicates the volume dependence of the temperature states ujl and the fluctuations 
Fl. The thermodynamic limit is computed in various ways, leading to the different phases 
of the Bose gas (normal, critical, and condensed). The difference between the critical 
and condensed phase is analysed by a special technique, namely by the interplay between 
the scaling of an external gauge breaking field, used to force the gas into an extremal 
state, and the speed at which the chemical potential converges. Using this technique, we 
obtain detailed information about the behaviour of the field and density fluctuations in the 
different phases of the Bose gas. As turns out, also the fluctuation distributions ( p..l|) and 
the scaling exponents in the fluctuation observables are very sensitive on the boundary 
conditions and on the way the thermodynamic limit is taken. We prove the existence 
of different regions in the space spanned by the scaling parameters of external field and 
chemical potential, in which the fluctuations are differently distributed and have a different 
degree of abnormality. The distribution functions also depend explicitly on the attractivity 
parameter of the boundaries. The distributions we obtain, are Gaussian or non-Gaussian, 
normal or abnormal depending on different choices of the scaling parameters. 

The main conclusion of this paper can be summarised as follows: details of the boundary 
conditions, the strength of the external field, and the thermodynamic limit have a vast 
impact not only on the condensation phenomenon, but also on the distribution functions 
of field and density fluctuations. These results show that the analysis of the thermodynamic 
limit should be done very carefully and its properties are of utmost importance for physical 
observable effects, which are detected at the level of the states as well as on the fluctuations. 

We remark that we take the thermodynamic limits of the states together with the volume 
scaling of the fluctuation observables. In this way, the dependence of the limiting distri- 
butions and the volume scaling exponents on boundary conditions can be made explicit. 
These limits are usually taken separately and boundary conditions are reintroduced with 
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special cut-off functions in the definition of tlie fluctuation operators 0] . This type of cal- 
culations was previously applied for the classical Curie- Weiss Model |T^, but as far as we 
know not yet for quantum systems. We obtain explicitly new distribution functions of fluc- 
tuations in our Boson model. In particular new non-Gaussian critical density fluctuations 
are derived. 



2 The Model 

2.1 The One-dimensional Eigenvalue Problem of the Free Lapla- 
cian 



We study the behaviour of a free Bose gas in a //-dimensional cube with edges of size L, 
centred around the origin A,^ = [— L/2, 

The first step in handling this multi-dimensional, many-body system consists of solving the 
basic one-dimensional one-body eigenvalue problem of the free Laplacian on the Hilbert 
space £^(Ai): 



[X) 



(2.1) 



with if) G C^(Ai), the two times continuously differentiable complex functions on Ai. We 
take the units ^ = 1 and consider a family of self-adjoint extensions of the Laplacian by 
restricting its domain using the following boundary conditions 



dip 
dx 



x) — aip{x) 



dtp / \ , / \ 
— [x) + aip{x) 

dx 



0, 




(2.2) 
(2.3) 



with cr G R, a parameter governing the elasticity of the boundaries. Implementing these 
conditions ( |2.1| )-( |27^ ), gives the spectrum and eigenfunctions. This is worked out in detail 
in Ijl^ , we mention here only the results. It can be proved that there is no continuous spec- 



trum and that there are infinitely many discrete eigenvalues (en)neN- The eigenfunctions 
are given by 



i^nix) = cos{^/e^x) if n is even, 
"Pnix) = sm{^/e^x) if n is odd. 



(2.4) 
(2.5) 



up to normalisation. The spectrum (e„)„ depends on the size of the box Ai, and on the 
elasticity parameter a. 

About the sign of the elasticity parameter, we distinguish the following situations: 
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Neumann boundary conditions (a = 0) 

The easiest case, is tlie case of Neumann boundary conditions i.e. a = 0. The system 
of equations (p.l|)-(p73D can be solved exactly and the spectrum is given by 

en=(^)', VnGN. (2.6) 

The wave function of the lowest energy level is a constant function. Physically, this 
means that the particles in the ground level are not attracted, nor repulsed by the 
boundaries, i.e. the situation cr = corresponds to perfect elastic boundaries. 

Attractive Boundary Conditions (a < 0) 

If 0" < 0, negative eigenvalues are present. Although an exact solution of the set 
of equations ( p.l|) -( p73D is impossible, graphical techniques can be used to deduce 
the following properties of the spectrum. If \a\L < 2 there is only one negative 
eigenvalue, and the spectrum satisfies the following useful spacing properties 

eo < -a^ < < ei < (^-j < ^2 < ( 1 < ^3 < ( J < £4 < • • • (2.7) 

If |cr|L > 2, the lowest positive eigenvalue becomes negative and there are two nega- 
tive eigenvalues. The spectrum behaves now as 

eo < -cT^ < ei < < < e2 < (^^^ < £3 < (^^^ < < ■ ■ ■ (2.8) 

Since we are especially interested in the thermodynamic limit with fixed elasticity 
parameter a, we always assume L\a\ > 2. The lowest eigenvalue is monotonically 
increasing to — a^, and the second eigenvalue is monotonically decreasing to — cr^, 
and these processes are exponentially fast: 

eo = -a^ - O(e-^l'^l), ei = -a^ + O(e-^l'^l). 



It follows from formulas ( |2.4| )-( p75| ) that the lowest eigenfunctions are given by cosh 
and sinh functions. Physically, this means that the particles in these orbits have a 
high probability to be close to the boundaries, or that the boundaries are attractive. 

Repulsive Boundary Conditions (a > 0) 

Using similar techniques as in the case of attractive boundary conditions, it can be 
deduced that there are no negative eigenvalues if cr > 0, and that the spectrum 
behaves as 

7r\2 /27r\^ /Stt^^ 



0<eo<(-j <ei<(^l <e2<(-^) < e-s < ■ ■ ■ (2.9) 
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The norm of the lowest eigenfunction decreases near the boundaries, or the bound- 
aries are repulsive. If cr t oo, we have Dirichlet boundary conditions, the eigenfunc- 
tions vanish at the boundaries and the spectrum can again be solved exactly, it is 
given by 

^" + ')"^', VneN. (2.10) 

These considerations about the one-dimensional case can be extended to more dimensional 
systems. Let us now consider the v dimensional Laplacians on Ky = [—L/2, L/2Y with 
similar domain restrictions, namely, replace (|2.2|)-(|2.3|) with 



dip 

— (x) = crip{x), X e dAy, 
on 

with ^ the inward normal derivative. The eigenvalues are now denoted by eLik), where the 
dependence of the eigenvalues on the size of the box is made explicit, k = {ki,k2, ■ ■ ■ ,ky) G 
N'", and the eigenvalues are given by 



eL{k) = J2^k. (2.11) 



2 = 1 



with the fcj-est eigenvalue of the one-dimensional free Laplacian on Ai (|2.6| )-( |27T0| ). The 



eigenfunctions ipk{x) G -^^(■^) ^-I's products of their one- dimensional components ( p.4|) - 



2.2 The Bosonic Many Body System 
2.2.1 The Hamiltonian on the Fock Space 

In order to be able to describe a gas of bosons in the box A = [— L/2, we use the 

standard techniques of second quantisation [15] and find the Hamiltonian on the Bose-Fock 
space 

Hl = Y^i^Lik) - fiL)a\k)a{k) - L^h (a(0)e-'<^ + a\<A)e"^) (2.12) 

k 

The index k runs through all vectors in N'', eL{k) is the energy associated with the level 



k (|2.11| ), < inffceL(/c) is the chemical potential, which determines the particle density 
in the system. The second term in ( |2.12[ ) is a external field term, breaking the gauge 
symmetry. It is added to recover one of the extremal translation invariant equilibrium 
states in the thermodynamic limit. This field determines the phase of the condensate, it 
scales with the volume with an exponent 7, which has to be chosen in the range 

-Z//2 < 7 < Z//2. (2.13) 
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If 7 < —'^I'^i its effect is to weak to cause a gauge breaking in the thermodynamic hmit, 
and if 7 > the field is too strong and causes an artificial gauge breaking, i.e. there 
is no non-zero critical density, /i is a positive constant. The creation operators a^(A;) and 
their adjoints, the annihilation operators a(fc) are defined by 

a(.)^/...^.(.)«(.). (2.14) 

where (i^{x) is the creation operator of a Boson at a; G Aj, and a[x) is the corresponding 
annihilation operator. The function tl}^^ G L'^i^Ky) is the eigenvector by the eigenvalue e2,(fc) 
( CT ). The local Hamiltonians are diagonalised in terms of quasi-particles, their 

creation and annihilation operators are denoted by 6^, resp. 6, and the relation between 
the annihilation operators of the quasi-particles and those of the bare particles is 

m = a(0)-^, (2 15) 

bik) = a{k), Vfc G N^\{0}. ^ ■ ^ 



This enables to rewrite the Hamiltonians ( 2.12| ) as 

eL(0) - fiL 



Hl = Y.^e,{k) - fiL)b\k)b{k) - j^j^^^^. (2.16) 

k 



We see that in terms of the quasi-particles the Hamiltonians Hl are essentially the free 
Hamiltonians plus an unimportant constant term. The volume dependence of the spectrum 
and the chemical potential are important in the remainder of this paper. 



2.2.2 Quasi Free Equilibrium States 

The finite volume equilibrium states or Gibbs states ul are now easily established using 
the expression in terms of the quasi-particles ( |2.16|) . The equilibrium states then coin- 
cide with the expression for the equilibrium states of the free Bose gas without external 
field [|I5|. These states are quasi- free states, the easiest way to characterise them is to 
use the truncated correlation functions u!l{- ■ ■)t, recursively defined, using the following 
expression, for all A,, i = 1,2, . . . creation or annihilation operators or combinations of 
them: 

uol{A, An)=J2Il • ■ • ' ^Jmh, Vn G N, (2.17) 

where the sum t E Vn runs over all ordered partitions r of a set of n elements in subsets 
J = {j(l), . . . , j(| J|)} G r. The truncated functions associated with the equilibrium states 
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ujl satisfy 



UJLibKf))T = 



^L(6n/i),K/2))T = $^/2(A:)/i(A;) 



1 (2.18) 



with /i, /2, . . . G LF'{Kv)i and 6" can be either 6 or 6^". Only the two-point functions are 
non-zero, the Fourier transform f{x) ^ f{k), k G N"^ is defined by 



/W = (^fc|/)= (2.19) 



3 Condensation 

In this section we prove the existence of different phases in the thermodynamic hmit 
of the system, i.e. we prove the existence of a Bose-Einstein condensate - a macroscopic 
occupation of the lowest energy level - when the chemical potential fi^ scales correctly. Here 
we take the thermodynamic limit {L f oo) with fixed chemical potentials fii and varying 
densities (as in |11]), rather than with a fixed density and varying chemical potential (as 



m 



12, TH]). But before discussing this in detail, let us first prove two lemmata. 



Lemma 3.1. Choosing the chemical potentials fii for each volume A,y = [—L/2, L/2Y 
equal to 

I^L = eL(0) - ^L--* (3.1) 



where the scaling exponent < a* = z//2 — 7 and po G IV~ , yield a non-zero occupation 
of the lowest energy level in the thermodynamic limit, i.e. with this choice for the series 
{fJ'L)L we have 

lim L-''ujL{a\0)a{0)) = po (3.2) 

and the gauge-invariance of the limiting state is broken, i.e. 

limL~''/^ujL{a{0)) = y^e't (3.3) 



Proof 

The density of particles in the lowest energy level in the equilibrium state ujl for a finite 
volume Aj, (|2.18|) is given by 

Po{L)=L-'uLia\0)a{0)). 
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This expression can be written in terms of the quasi-particles using the relations ( p.l5| ) 
and the expression for fx^ ( |3.1|) , i.e. use a(0) = 6(0) + L'^/^y/poe"^ to find that 

PoiL) = L-'LULib\0)biO)) + L-'/'^^o K(a(0))e-''^ + ^L(at(0))e'*) - po- 

The first term is of order 0{L~'^^^~"^), this can be seen using the two-point functions ( p.l8|) 
as follows 

^-u 1-0{L-') 

+ +0(L-2"*) 

= 0{L''*-'') = OiL'"/^-^). 

And by ( |2.13| ), this term vanishes in the thermodynamic limit. We use this observation to 
establish a bound on the difference between po{L) and po 

< po{L) - L-'/^^oiM^^me-"^ + ujLia\0))e"^) + po < 0(1-"/^-''). 

Using the Schwarz inequality L^''/'^\LUL{a{0))\ < L^'^iUL{a'< {0)a{0)) = a/ PoiL) the above 
expression is transformed into 

Po{L)-2^o^/M^ + Po < 0(1-"/'-^) 

and hence, 

\^/p'o-VML)\ < 0(L-/^"^/2) 

l(v^- v/p^)(v^+v/M^)l < 0(L-/^-^/2) 

\po-Po{L)\ < 0{L~'^"~^l^). 

Indicating that po{L) converges to po in the thermodynamic limit. 

Also the second part of the lemma ( |3.3[ ) can be proved by means of the Schwarz inequality, 

L-n^L(&(0))|' < L-''ujL{h\m^)) = 0{L~^/'-^). 
Using the relations (|2.15|) and the expression for the chemical potential, this yields 

hm L-'/'uMO)) = v/p^e'<^. (3.4) 

L— >oo 

□ 
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The result ( p. 21 ) is not yet a complete proof of Bose- Einstein condensation, it should also be 
proved that the total particle density is finite with this choice for the chemical potential. 
Only then one can speak of a macroscopic fraction of particles condensed in the lowest 
energy level. 

The total density py{L) in a volume in the equilibrium state uj^ is easily derived from 
the relations (|2.15|) and the two-point functions (|2.18|) 



k k^O ^ 

where zl is the activity or fugacity zl = e^^^, and the sum over k runs over all k G N'^yjO}. 
We formulate now the well-know result 

Lemma 3.2. Take zl z < lim^, e^^^ then 

-J^ = A- V(z//2, z) (3.6) 

where A = tt\/J3 is the thermal wavelength. The function 

J{v/2,z) = [ d^x 



ze 



-J. 2 

-.^i^u — ze 

is the Jonquiere function. This result is valid for any type of boundary conditions specified 



Proof 

Based on the spacing properties of the eigenvalues (|2.6|) - (|2?T0|) and the convergence of 
Riemann sums to Riemann integrals. A proof can be found in 



□ 

The Jonquiere function J(z//2, z) can be expressed as 



oo 

Z 



^('^A^) = Ew^ (3-7) 



n=l 



It is an analytic function of z in the cut-plane, the cut being from 1 to oo along the positive 
real z-axis. In the limiting case 2; — > 1, it converges to the Riemann zeta function ({i>/2), 
this is finite for dimensions higher than three or u > 3 ||16|| . 

Combining the results of the two lemmata we can analyse the different phases of the Bose 
gas. If we take the thermodynamic limit such that the series {pl)l converges to a certain 
value yLi, strictly lower than the limit of the lowest eigenvalue n < ji^. = lim^^oo ^^(0), there 
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is no macroscopic occupation of the lowest energy level (cfr. lemma and the total 



density converges to a finite value (lemma p.2|) . This situation is called the normal phase. 

The critical and condensed phase can be reached if we take a series converging 
to the maximal value /i,,,, provided that the total density ( |3.5|) is finite in this limit, i.e. 
in dimensions i/ > 3 if we have a > boundary conditions (where = 0), and in al 
dimensions for o" < boundary conditions (where = —ua"^ < 0). As said before, we 
study in this paper the critical and condensed phase in the thermodynamic limit {L | oo) 
taking series of the form 

/iL = eL(0)-cL-", (3.8) 

where c > some positive constant and a > a scaling exponent. Clearly, these series 
( p.8|) converge to fi^ = hm/, 6^(0). Depending on the scaling exponent a, we have a different 
phase in the thermodynamic limit. 

If the convergence is slow, i.e. if the scaling exponent a is chosen in the range 

< a < a, = Z//2 -7, (3.9) 

we end up in the critical phase. The total density converges to its critical value A^*^ J(i//2, 2;^,), 
with = e^^* . There is no macroscopic occupation of the lowest level, since the conver- 
gence /iL — > /i^, is too slow (lemma 



If the convergence is sufficiently fast, i.e. if the scaling exponent a (|37^ ) is large enough, 

a = a, = z//2-7, (3.10) 



i.e. if /i/, converges to /i* as specified in lemma p.l| , the total density can reach arbitrary 
values above the critical density Jiv j1^z^\ There is a non-zero density of the con- 
densate po- This condensate density is determined by the constants c ( |3.8| ) and h of the 
external field (|2.12|) , i.e. po = ^ jf?-, cfr. equation (|3.1|) . 

The situation with \il = — cL'^" with a > z//2 — 7, is not physically meaningful 
since we end up in a situation where the density of the condensate Pq{L) diverges in the 
thermodynamic limit. 

Let us summarise this in the following. 

Theorem 3.3 (Phases of the Free Bose Gas). Fix a temperature (3 , the thermodynamic 
limit of the (3 -equilibrium states ujl, i.e. the limit {L f 00) with (/i^ fi), exists if the total 
density 3i) converges to a finite value p, and 

• p < Jiy l2, z^) , if pl ^ p < p^ = limiei(O) (Normal Phase), 

• p = J{u /2, z^), if pi = £^(0) — cL^"" with c > and < a < i^/2 — 7, and if 
J{i'/2,z) converges (Critical Phase), 
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• p = pq + \ " J{u/2,z^), if Pl = eL(0) — -;^L °' with = z//2 — 7, provided that 
J{v/2^z) converges (Condensed Phase). 

For attractive boundary conditions, the condensation phenomenon takes place in any di- 
mension, but due to the special form of the wavefunctions of the lowest energy levels, 
the condensate is localised near the boundaries, and the condensation is a pure surface 
effect. For o" > 0, the condensation is a bulk phenomenon and only takes place in di- 
mensions higher than three (z/ > 3). More detailed discussions of these different types of 



condensation can be found in |[TT|,p!2|,p!5[] 



It is clear that the strength parameter of the external field 7 plays an important role in the 
condensed and critical phases. If 7 is small, the a has to be larger, or p^ has to converge 
faster to zero, in order to have condensation; for larger values of 7, condensation is already 
present at a slower convergence rate ioi pi. A value 7 = z//2 provokes an artificial gauge 
breaking, and there is always condensation at any density. If 7 < — the nature of the 
phase transition changes, the condensate is no longer of a well defined phase, and we do not 
find a single extremal state in the thermodynamic limit, but a mixture. Hence clustering 
is absent and we can no longer analyse the fluctuations. The interplay between 7 and a 
plays also its role below where we analyse the field and density fluctuations. 



4 Field Fluctuations 

In this section we study the scaling behaviour of field fluctuations, or fluctuations of oper- 
ators of the form 

At = l={a{k) + a\k))- (4.1) 
A- = ^{a{k)-a\k)), (4.2) 

with k G N*^, different modes. The local k-mode field fluctuations are defined by 

FlA^^) = {At - MA^)) (4.3) 

5 is the scaling exponent and should be chosen such that the limiting characteristic function 

^(F,(A|)) : t ^ ifiFsiA^m ^ hm c.z.(exp(itFi,,(A^J)), (4.4) 

is non-trivial [0. The limits we consider here are limits of series of local fluctuations 
FL,5iAt^), where the vectors fc^ G N*^ are chosen such that we can associate with the series 
(k^)^ a series of eigenvalues {eL{ki))j^, converging to a certain value E G {— z/cr^} U [— (i/ — 
1)0"^, 00) in the limit L ^ 00 . The scaling exponent 6 depends on E and on the phase of 
the Bose gas. 

In the case of field fluctuations, the limiting characteristic function can explicitly be ob- 
tained, due to the quasi-free character of the states ul- 
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Lemma 4.1. Field fluctuations are Gaussian 



ipiFsiAi))(t) = lim exp ( ~L-'' (l + 2uj L{b^ {kL)b{kL)))] • 

L^OO \ 4: J 



(4.5) 



Proof 

The characteristic functions ( ^.4| ) can be written in the following expansion 



°° (it)"' 

(exp(itQ)) = exp — -ujl{Q, <5, • • • , Q)t 



(4.6) 



n=l 



n times 



with ljl{Q, Q, . . . , Q)t, the n-point truncated correlation functions ( p. 171) . In the case of 
field fiuctuations, i.e. if we substitute Fi^si^t) Q (|l-6[), due to the properties of the 
equilibrium states ul (|2.18|) , only the second order term in this expansion is different from 
zero. This term can be rewritten in terms of the quasi-particles as 



J -26 

= -—{l + 2ujL{h\kL)h{kL))) 



and this yield the explicit form of the characteristic function ([4.5| 



□ 



In spite of the fact that field fiuctuations are Gaussian, they can be abnormal, in the sense 
that there has to be a non-zero scaling exponent 5, in order to have non-trivial fiuctuations. 

Theorem 4.2 (Field Fluctuations). The limiting characteristic functions of the k-mode 
field fluctuations 

^{F{A%){t) = hm (exp(itF^,5(A± ))) (4.7) 
with eLikh) E, tend to non-trivial distributions if 

• 6 = in the normal phase [fiL — > < /i* = —ua'^) 

^{FiA^m = exp(-^ coth(/5/2(i? - /i))) (4.8) 

• 5 = in the critical and condensed phase if E ^ — z/a^ 

y.(F(A|))(t) =exp(-^coth(/?/2(E + z/a2)) (4.9) 
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• 5 = a/2 in the critical and condensed phase {^l = — cL for all a : < a < 
Z//2 — 7 and if tLlk^) E = —ua"^ 

^{F{A%^Mt) = exp{-^^J. (4.10) 

In the condensed phase, i.e. if a = v/2 — 7 this expression can also he expressed in 
function of the condensation density po, i.e. 

<^(^«.0)W = exp(-^^). (4.11) 

Proof 

The two-point function appearing in the expression for the distribution of field fiuctuations 
( |4.5|) is given by 

MbHkM = z^"'",^';^''}':"^" V (4.12) 

1 - exp{-(j{eL{kL) - PL)) 
This function converges in the hmit {L | 00) if 

hm \eL{kL) - Pl\ > 0. 

L— >oo 

This situation occurs in the normal phase for all eilki) — ^ E, but in the critical and 
condensed phase only for 62,(^2,) E ^ —va^. It yields normal fluctuations, i.e. (5 = 0, 
and the explicit form of the distributions ( |4.8| ) and ( ^.9| ). 

Let us now consider the case E = —ua^ in the critical or condensed phase. If E = —ua"^, the 
vectors ki should be in {0, 1}*^ for sufliciently large L, only in that case there is convergence 
^liki) —lycr"^, and this convergence is exponentially fast. In the critical phase pi also 
converges to — z/cr^ (^.81), and the expectation value appearing in expression ( |4.12| ) diverges. 
An extra scaling L"^ is necessary in order to obtain flnite variances. The highest order 
term of ( [4.12|) behaves as 

p-f3{<^L{kL)-IJ-L) 

2L-V(6t(fc,)Kfc.)) = 2L--^--^^,^^^ 

l-/?cL-" + 0(L-2") + ... 



l-l + /?cL-° + 0(L-2") + ... 
2L-2^(^L" + 0(1) + .. 



Hence, a scaling 5 = a/2 is needed in order to avoid divergences or trivial distributions. 
The distribution then converges to 

¥'(i^«..))W = exp(-^^). 

In the condensed phase, i.e. if a = z//2 — 7, we can replace c by an expression depending 
on the condensate density and the external fleld h: ^ = ( [4.111 ). 

□ 
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5 Density Fluctuations 



The local density fluctuations -Fl,<5(^o) are defined by 

FL,5iNo) = L-"/'^-^ [ dx a^{x)a{x) - coiia^ {x)a{x)), (5.1) 

J A 

or equivalently 

Fl^No) = L--'^-' J2 ^^(pHp) - MaHp)a{p)). (5.2) 

it turns out these density fluctuations can be divergent in the critical and condensed phases, 
and suitable non-zero scaling exponents 6 have to be added. 

Instead of using such an extra scaling factor, one can introduce modulated fluctuations. 
Mostly one chooses a cosine function as modulation function, and one introduces the so 
called fc-mode fluctuations by: 

Fk{A) = L-^l^ ! dx {tM) - ^L{rM))) cos{kx). 



Inspired by P], where the Goldstone phenomenon was studied in interacting Bose gases, 
we consider here fc-mode density fluctuations given by 



FL^Nk) = L-'^l^-']^ Yl {^^(pHp + k) + a^ip+ k)a{p) 



UL{a\p)a{p + k) + a\p + k)a{p))) . (5.3) 



2 

p 



These observables can be considered as exchange correlations between particles of mo- 
mentum differences k, but also as modulated density fluctuations. In the case of periodic 
boundary conditions, which was considered in |P, this form of modulation coincides with 
the fc-mode cosine-modulated fluctuations. Here, in the case of elastic boundary conditions, 
these fluctuations (|5.3| ) correspond to a more complicated form of modulation 



FL,s{Nk) = L-"/^-' J^^^jj^y aHx)a(y) (^^^p(x)^p+fc(y)/2^ + h.c. 

The modulation function ■ ■ ■ is here no longer of a delta-function type, i.e. 6{x — 
y) cos{kx), but it is a more smeared out function. 

These fluctuations (prT|)-(pl3|) have no longer a priori a Gaussian distribution. Nevertheless, 
we are still able to compute rigorously and explicitly their distribution functions as a 
function of temperature, density, external field strength and boundary conditions. 

For the sake of compactness of the paper, we will not discuss in this section the case where 
k depends on the volume k = ki (as in Theorem [4.2|) , but take k constant. The more 
general situation however is easily obtained from this case. 
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5.1 Quasi-particles density fluctuations 



The distributions of the density fluctuations of the bare particles is related to the distri- 
butions of similar density fluctuations in terms of the quasi-particles ( p.l5 ). The first step 



in calculating the distribution of density fluctuations consists in finding the distributions 
of the density fluctuations of the quasi-particles, 

Fl,5{K) = L""!^-']^ {bHp)KP + k) + b\p + k)b{p) 

p 

- ooL{bKp)b{p + k) + b\p + k)b{p))], (5.4) 



where the accent in N'^, is added to make a distinction with the particle density of the bare 
particles. 

5.1.1 A; 7^ quasi-particle density fluctuation 

First we consider the A; 7^ density fluctuations. In this case, the expectation value 
appearing in ( |5.4| ) is zero and may be left out. If we want to calculate the characteristic 
function 

ifit) : t ^ lim UlU'^^'^^^'A 

L^co \ J 

we can use expansion ( [4.6|) in the n-point truncated correlation functions: 

These n-point truncated functions can be calculated, but first we need to know something 
about the expectation values of the powers of F{Nk). 

The expectation value of powers of the quasi particle density fluctuations. Con- 
sider the function uJilF{N'j.Y) with F{Nl) as in (p^). This function can be expanded in 
2" functions, if every factor F{N'j.) is splitted into two parts 



L-''/y2Y,bHp)b{p + k)j + (^L-'/y2^b^{p + k)b{p)j . 



(5.5) 



p 

Those 2" terms are all of the form 



J2 u;L{b\pi)b{p, + k)...b\p^ + k)b{pr,)) (5.6) 



2 

Pl,-,Pn 



where n factors appear being either b^{pi)b{pi + k) or b''{pi + k)b{pi), Vi = 1, 2, . . . , n. Such 
a term (|5.6|) can be represented as a configuration of n symbols on a circle, e.g. see (Fig. ||) 
with n = 6. 
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Figure 1: A possible configuration with 6 sites and 3 cycles. 



We draw a circle and add n points, choose a starting point and count the different sites 
clockwise. To every site 1, 2, 3, . . . , n we add a symbol being either o or x. We draw a 'o' 
at site i, for a factor h'^{j)i)h{pi + /c), and a 'x' for a factor h\pi + k)h{pi). In the example 
(Fig. |1|), we represent the monome 

ujl {h\pi + k)b{pi)b\p2 + k)b{p2)b\ps)b{p3 + k) 

Pi,... ,P6 

b\p4)b{p^ + k)b\p5 + k)b{p5)b\pQ)b{pe + k)) . (5.7) 

Using the quasi- freeness of the states ( |2.18| ), these functions ( |^ ) can be expressed as a 
sum of terms consisting of n products of 2-point functions only. All partitions into ordered 
pairs appear once in this expansion. A possible term in this expansion of products of 
two-point functions is now represented by a directed graph on this circle, connecting the 
different sites with each other, obeying the rule that on every site there must start an 
arrow and must end one. Every arrow corresponds to a two-point function. Such a two- 
point function is constructed by combining the creation operator from the site where the 
arrow starts with the annihilation operator from the end-point. The order of the operators 
is defined by the order of the sites 1,2, ... ,n. A loop is also possible, then we combine 
the creation operator with the annihilation operator from the same site, e.g. in the above 
example (Fig. HD, there is a loop at site 6 , this yields a factor u!L(b^ {pe)b{pQ + fc)). Since 
k ^ 0, such loop-factors are zero by gauge-invariance of the states ul ( p.l8|) . A graph on a 
circle consists of separated connected graphs on subsets or cycles, e.g. the graph in (Fig. ||) 
consists of a two-point cycle, a three-point cycle and a one-point cycle or loop. The total 
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term represented in (Fig. |1|) reads 

UJL{h\pi + k)h{p^)) UJL{b{p3 + k)b\p5 + k)) UJL{b{pi)b\p3)) 

pi,--- ,P6 

uJL{b\p2 + k)b{p, + k)) UL{b{p2)b\p,)) UJL{b^{p,)b{p, + k)) . (5.8) 

A number of useful properties can immediately be deduced from this representation: 

• Every (non-zero) cycle has only one summation index. 

Every two-point function yields a relation between the summation indices of the 
creation and annihilation operator. All types of two point functions appearing in 
these expressions (|5.8| ) are zero unless the indices of the creation and the annihilation 
operators are the same. This yields a linear relation between the summation indices 
Pi,P2, . . . in the cycle, and hence only one summation is free. 

• A cycle containing not the same numbers 'o' as 'x' are zero, consequently cycles over 
an odd number of sites are always zero. 

The sum of the indices of the creation operators minus the sum of the indices of 
the annihilation operators in a cycle must be zero, otherwise there will always be a 
two-point function where the indices of the two operators are different, and such a 
factor is zero. In cycles where there are not as many symbols 'o' as 'x', this sum is 
always different from zero. 

Truncated functions. The following step in the calculation of the characteristic func- 
tion, consists in the calculation of the ra-point truncated correlation function itself. 

Lemma 5.1. All odd truncated functions vanish. 

UL{FiNi),...,FiNi))^ = 0, VneN. 
^ ^ ' ^ 

271+1 factors 

Proof 

First note that the one-point function vanishes, 

u;L{F{N',))T = L-''/y2Y,uJL{b\p)b{p + k))+u;L{b\p + k)b{p))=0 

p 

due to the gauge invariance of the states ul- 

Consider now the 2n + 1-point truncated function, and suppose that all m-point truncated 
functions where m is an odd number less than 2n + 1, vanish. Using the definition of the 
truncated functions (|2.17|) , the 2n + 1-point truncated function is written as 

(5.9) 
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where the sum in the second term on the rhs runs through all partitions r eV in two or 
more subsets J of a string of 2n + 1 elements. Each term in this sum on the rhs contains 
at least one factor with a truncated function over an odd number of points m < 2n + 1, 
and such a factor is zero by induction hypothesis. 

Also the first term ujl{F{N'i^)'^^^^) vanishes. Expand this term using the definition of 
F{N'i^) ( p.4| ), as in ( |5.6|) . All possible configurations in terms of two-point functions (Fig. 
|l]), contain at least one cycle with not as many symbols 'o' or 'x', because the total number 
of symbols is odd, and such cycles vanish. Hence there are no configurations which are 
non-zero and u l{F (N^.)'^"''^^) = 0. Hence, also the 2n-|-l-point truncated function vanishes, 
and by induction, all odd truncated functions vanish. 

□ 

Lemma 5.2. All even truncated functions can be written as the sum over all configurations 
with two-point functions (as in Fig. \^ containing only one cycle connecting all sites. 

Proof 

Consider first the 2-point truncated function 

ujmK),F{N',))^ = u;,{F{N',r) -uJriF{N',))u;,{F{m 

= L-'^i ^L{b{pi)b\p2)) u;L{b\pi + k)b{p2 + k)) 

P1.P2 

+ ujL{b{pi + k)b\p2 + k)) UL{b\p,)b{p2)) (5.10) 

Using the diagrammatic representation of these functions (as in Fig. |I|), we see that 
uJilF{N'f,), F{NQ)rp can be written as the sum over all (non-zero) diagrams with two sites 
containing only one cycle, cfr. (Fig. ^). 




Figure 2: The non-zero diagrams for ujL{F{N'f.Y) 

Consider now the 2n-point truncated function 

u,{F{K\ . . . , F{N',))t = MFiKrn ^l{F{K), . . . , F{K))r, (5.11) 
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and suppose that all 2m-point functions, with m < n are of the prescribed form. Further- 
more all odd truncated functions vanish (lemma |5.1| ), hence the sum in the second term is 
a sum over products of diagrams with one cycle and an even number of sites. 

The first term uo ilF {N'/^)^'") can be rewritten in terms of diagrams (as in Fig. 0). We use the 
fact that diagrams containing several cycles can be written as the product of diagrams over 
less points containing only one cycle, because all cycles are independent. Take a certain 
partition r G P of the 2n sites into subsets of | Ji|, | J2I, . . . ,|^r,„„^| sites, and consider 
now all possible arrow- diagrams where there is for every subset Ji, J2, . . . G r just one 
cycle, connecting all the points in that subset. By distributivity we can rewrite this as the 
product over all subsets Jj of the sum of all possible configurations with one cycle on a 
diagram with | Jj| elements. 

Since all cycles on diagrams with an odd number of points vanish, we only have to deal 
with the partitions r in subsets all containing an even number of points. In the case of a 
partition r in two or more subsets, we can apply the induction hypothesis, i.e. suppose that 
the sum over all possible diagrams with 2m < 2n points with only one cycle is equal to 
the 2m-point truncated function. This yields that the term corresponding to the partition 
T can be written as the following product 

l[ujLiFiK),...,FiN',))^. (5.12) 
Hence, the 2?T,-point truncated functions ( |5.11|) can be written as 

^,(F(ivo,...,F(iv^))^ = +En^^(^™'---'^™)^ 

diagrams with one 2n-cycle T" 

v 

diagrams with one 2n-cycle 

Also the 2n-point truncated function is of the prescribed form and by induction, this is 
valid for all even truncated functions. 

□ 

Now we are ready to calculate the distributions of the quasi-particle density fluctuations. 
Theorem 5.3. The k ^ quasi-particle density fluctuations are 

• Gaussian and normal 

v{F{Km=eM-t\{z)/^) (5.13) 
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with 

in the normal phase, and in the critical and condensed phases (fi^ = — cL~°' , 
cfr. (^)) for k ^ {0, l}^ or for k e {0, 1}^\{0} and a < z//2. 

• Non- Gaussian and normal 

^(FiN'Mt) = e-*^«^)/^ ( i— ^) (5.15) 

VI + (4) / 

with C(2*) as in §J^, = e"^^'^", and a{k) = 2^''-^^\ for k e {0,1Y\{0} and 
a = u/2. 

• Non- Gaussian and abnormal (5 = a — v /2) 

^imm) = ( — ^] , (5.16) 

with a{k) = 2(^-'='), for k G {0, 1Y\{0} and a > z//2. 
Proof 

Applying the expansion ( [4 .61) to the functions F{Nl.) yields the following expression for 
the characteristic function 

u;,{e^^^(K)^ = exp J2 ^^^(^(^^)' • • • ' ^(^^))t (5-17) 

n>l 

where ul{F{N'}^), . . . , F(A^^))^ are the n-point truncated functions. 

Normal Phase: Lemma and lemma ^]2| about the truncated functions learn that all 
odd truncated functions vanish, and that all 2n-point truncated functions, 2n > 2 are 
vanishing, this last fact is easily seen as follows. From lemma f).2[ the 2?2-point truncated 
functions could be written as a the sum over all possible cycles with 2n points. Such a 
terms are of the form 

J2 ^L{h\p)h{p)) UL{h{p + 2k)h\p + 2k)) ■ ■ ■ (5.18) 
p 

with only one summation index, and 2n two-point functions, ujL{h{p + jk)b^(p + jk)), j = 
0, 1,2, . . . or with the operators in different order. All these factors are bounded in the 
normal phase, and using the techniques of lemma |3]^ it is easy to see that the sum 

Yl ^L{h\p)h{p)) ujL{h{p + k)h\p + k))--- UL{h\p + ?>k)h{p + 3A;)) 



20 



converges to a finite integral. But we have an extra scaling factor L ^'^ ^ which makes 
the 2n-point truncated functions, 2?2 > 2 behave like ©(L^*^"^^^/*^). 

Hence, the only extensive term is the contribution of the two-point functions ( |5.1CI| ): 
lim ujlFiN'^), F{N',))^ = lim L"'^^ uj^iKp^pip^)) + k)b{p2 + k)) 

Pl,P2 

+ ujL{b{pi + k)h\p2 + A;)) ujL{h\pi)h{p2)) 




= az)/2, (5.19) 

with < 2; < e"^'^'^^ for the normal phase and A = T^\f^i the thermal wavelength. 

One can check that there exists a constant M(z) > such that the 2n-point truncated 
functions, are bounded by 

2r)l 

■ (2n- 1)! L-'^("-i)M(2)", (5.20) 



n\n\ 
and 

lim y M!!M(2n - 1)! L-^^^-"-^ MizY 
L^oo^ 2n! n\nV ' ^ ' 

n>2 

converges to zero as 0{L~^). Hence, the error series can be controlled and the distribution 
function converges to (|5.13|) . 

Critical and Condensed Phase: In these phases, the contributions of terms like u L{h{q)h'^ {q)) 
diverge if g G {0, 1}" . The rate of divergence depends on the scaling exponent (|3^), we 
have 

^L{h\q)b{q)) = j^L^ + 0(1) = UL{h{q)h\q)) , Vg G {0, 1}^ (5.21) 

A 2n-point truncated function consists of terms of the form ( ^.18|) , with one summation 
index p. Only the first terms in this sum over p G N, with p G {0, 1}^ , in such a term 
contain diverging factors ( p.21| ). 



If A; ^ {0,1}'^, the only diverging factors are of the form oj L{h{j))l)^ {pj) or uj {j))h{j>)) , 
p G {0, 1}" . All factors with indices p + jk, with j = 1,2, . . . and arbitrary p, are finite. 

A term represented by a cycle which connects alternating symbols x and o, contains n fac- 
tors of the form co'j;^( 6 (p)6"l'(p)) or uiip^ {p)b{p)^ and n factors of the form ti;i(6(p -|- k)b''{p + /c)) 
or LJi{b^{p + k)h{p + fc)) . In cycles with arrows connecting two symbols of the same kind x 
or o, some of these factors are replaced by factors ujilpij) + jk)h'^{p + j/c)) with j = 2, 3, . . . 
depending on the number and the order of such an arrows. Hence, the maximal number of 
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diverging factors in a term of the 2n-point truncated function is n if /c ^ {0, 1}^, and the 
terms containing the largest diverging factors are of order 

Since 0<a<z//2 — 7<z/, these terms are vanishing, and the fc-mode quasi-particle 
density fluctuations, with k ^ {0, 1}'' are normal (no extra scaling exponent 5 7^ needed) 
and Gaussian as in the normal phase 



lim ulU'^^'^'A = e-*'«-)/^ 



with C(-2*) as in ( |5.14| ) and z^, = e 



The situation is different for the k-mode fluctuations, with k G {0, l}'^. Now not only 
factors with index p are diverging but there are also diverging factors with index p + k & 
{0, lY- The factors with indices p + jk, with j = 2, 3, 4, . . . which appear in cycles where 
arrows between symbols of the same type are present, are bounded since jk G {0,jY, and 
thus p + jk ^ {0, lY- The diagrams with the highest number of diverging factors are those 
where alternating symbols x and o are connected. They have then 2n diverging factors of 
order 0(L°). Such terms become extensive if a > v/2. For a 2n-point truncated function 
there are (2n)!/(n!)^ diagrams with an equal number of both symbols. For such diagrams 
there are n\{n — 1)! possible cycles which connect all points alternating the two types of 
symbols. Hence there are (2n)!/n terms of leading order. In leading order they are all 
equal to 



L-nu-2nS ^ 2n /l\2n 

p,p+ke{o,i}'' 

The sum J2pp+k€{oiy yields only a numerical factor a{k) = 2^"^^^^ 
Hence, inserting this in expression ( p.l7| ) yields 

lim .;,(e'*^(^^)) = lim exp ^^.(F(iVO, . . . , F{N',)), 

\n>l 

= -p(^-(^)E^v^2^) j 

= exp I (t(A;) V- (^— ^ 
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In the case a = z//2, these terms are extensive as well as the extensive term from the case 
a < Z//2 and the distribution is equal to the product of the distribution in the two other 
cases. 



lim cuife'" ^^'fe 

L— >oo V 



with C(-2*) as in (|5.14|) . This distribution is normal {6 = 0), but non-Gaussian. Finally, 
remark that also in the critical and condensed phases the series over the subdominant 
terms can be controlled using similar arguments as in (|5.20|) . 



□ 



5.1.2 Unmodulated quasi-particle density fluctuations 

The situation is similar for A; = quasi-particle fluctuations: 



F(iV^) = L-"/' Yl b^(P^b(P^ - ^L{h\p)h{p)) (5.22) 

pew 



In order to calculate the characteristic functions, we take again the expansion in terms of 
truncated correlation functions (^]6|), 

e'*^W) = exp ^c.,(F(iV^), . . . , F{K))^ , 

n=l 

and we look for an expression of the n-point correlation functions. 

First we need to understand the form of the expectation value of powers of F{Nq). By 
the quasi-freeness and the gauge invariance of the states, the monome ujl{F{N'qY) can be 
decomposed into a sum over all pair-partitions, where every pair-partition corresponds to 
a term consisting of a product of two-point functions. Just as for the /c-mode fluctuations, 
we can visualise this in a diagrammatic representation of ujl{F{Nq)^), (e.g. Fig. ^). 

Again we draw n sites on a circle, and label them from 1 to n. Each site corresponds to 
a factor L~'^/^ l)^{pi)h{pi) — ujL{b'^ iP'i)b{pi)y Each pair-partition can be represented by 
a directed graph (Fig. ^). In every site starts an arrow and ends an arrow, every arrow 
represents a two-point function, this two-point function is constructed in the following 
way, take the creation operator from the starting point of the arrow, and combine it with 
the annihilation operator from the endpoint, the order of the operators in the two-point 
function is imposed by the order of the sites, the operator with the lowest site number 
comes first. Loops are now not permitted, the effect of the subtractions of the expectation 
values uji,(b^ {pi)b{pi)) in the definition of the fluctuation observables (|5.22|) , lies just in 
the cancellation of the terms or graphs containing loops, i.e. factors to L{b^ {pi)b{pi)) . Every 
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Figure 3: A 6-point diagram with a directed graph consisting of a 2-point cycle and a 
4-point cycle. 

graph consists of a bunch of independent connected subsets or cycles. Since every two-point 
function is only different from zero if the indices of both operators are equal, all summation 
indices of the points in a cycle are equal, and there is only one effective summation index 
for every cycle. 

Lemma 5.4. n-point truncated functions uJilF (Nq) , . . . , F{Nq))j., with n > I, can be writ- 
ten as the sum over all diagrams over n points with one cycle connecting all n points in 
the diagram. 

First note that the one-point truncated function is zero, 

ul{F{K))^ = u:iF{N'^)) = 0. 

The two point function reads 

uJiiF{N',),F{N'^)^ = ujL{F{N',f) -u;iiF{N',))ujL{F{K)) 

Pl,P2 

This is the term represented by the only diagram over two points with one cycle connecting 
both points. 

The rest of the proof uses induction. Consider the n-point truncated function and suppose 
that all m-point truncated functions, with m < n are of the prescribed form. The n-point 



24 



truncated function is defined as 

u;i(F(iV^), . . . , F(iV^))^ = ^L(F(iV^)") -EH ^^(^(^o), • • • , mo))T , (5-23) 

rev Ji^T 

where the sum over V is over all partitions r into two or more ordered subsets Ji, J2, ■ ■ ■ G r. 
Since the one-point truncated function is zero, partitions containing singletons may also 
be omitted. The first term is the expectation value of the n}'^ power of F{Nq). It can be 
written in terms of diagrams. Using similar arguments as in lemma |5.2| and the induction 
hypothesis, one can write it as 

^mKr) = E + E n ^^(^(^0), • • • , mo))T , 

n-point diagrams with one cycle tSP' JiSr 



as the second term in this equation is equal to the second term in ( ^.23|) , we conclude that 



u;i(F(iV^),... ,F(iV^))^= E 

71-point diagrams with one cycle 

or the n-point truncated function can also be written as a sum over all possible diagrams 
with n-points connected through one cycle, and by induction all truncated functions are 
of this form. 

□ 



Theorem 5.5. k = quasi-particle density fluctuations ( [5. 23^ ) are 
• Gaussian and normal 



^(F(iV^))(t)=exp(-t2C(^)/2) (5.24) 



with ({z) as in \5.1^ ),in the normal phase and in the critical or condensed phases 



(with fiL = eL(0) — cL if a < vjl; 
Non- Gaussian and normal 



,(FW))W^e-«.../^(^-^j (5,25) 



with C(-2*) CLS in (|5. 1^, in the critical phase if a = vjl; 



Non- Gaussian and abnormal (5 = a — v /2) 



in the critical or condensed phase if a > v/2. 
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Proof 

In lemma ^.41 we learned that an n-point truncated function can be written as a sum of 



terms of the following form 

peN" 

where the n factors are either h'^{p)h{p) or h'^{p)h{p). In total, there are (n — 1)! different 
cycles connecting all points in a diagram of n points, hence there are (n — 1)! terms of this 
form. We can analyse them as follows 

/2 u:L{h\p)h{p))... + L-^l' J2 ujL{h\p)h{p))--- (5.27) 

In the normal phase all factors are bounded, hence since the first term is a finite sum over 
bounded terms, due to the scaling factor L""''/^, it vanishes in the thermodynamic limit. 



Using the results of lemma 3.2 



lim L-"" V ujL{h\p)h{p)) = X-" dx - 

L — >CxD ' ^ I 1 

P^{0,1}- 



ze 



ze 



and since all two-point factors are bounded, also the second term in ( |5.27| ) is finite. The 
extra scaling factors L~"^/^ make that these terms vanish if > 2. Hence only the two- 
point truncated function contains an extensive term, it reads 



limL-^ J2 ujL{bHp)b{p))^L{b{pp{p)) = X'^ [dx . , ^_ 



ze \ 2 ze 



p^{0,l}- 



with ({z) as in ( |5.14]) . 

In the critical and condensed phases, the analysis of this second part of ( ^.271) is analogous, 
but now, the first part contains diverging factors, cfr. ( p.2H ). If a < z//2 the scaling 



exponent L ^1"^ dominates and the first term of (|5.27|) vanishes. The situation is thus 



similar to the normal phase, we have normal fluctuations (5 = 0) and the only dominating 
term is the two-point truncated function, 

^im L-'^ ^ ^L{h\p)h{p)) uj{h{p)h\p)) = ((z,), with z, = e"^'^'^' < 1. 

If a > Z//2, the first part in ( |5.27| ) becomes dominant, the fluctuations become abnormal 
and non-Gaussian, a scaling factor {S = v jl — oi) is necessary in order to obtain non-trivial 
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distributions. Using the relations ( |5.21| ), the expression for ( [5.27| ) is in leading order equal 
to 

J2 a;46t(p)5(p))... = 2"-L- + 0(L™). 
This yields the following expression for the characteristic function of F{Nq) 



ewj:'-^u.iiF{N',),...,F{K))^ = exp5:i-f(n-l)!2^ 

n=l ' n=2 ' ^ 



oo 

exp (2^ 



-it/l3c 



n=2 

2" 



it 

l3cj n 



1 - it /13 c 



As in the case for 7^ quasi-particle fluctuations (|5.20 ), it can be shown for A; = 0, that 
the series over the subdominant contributions vanish in the infinite volume limit. 

If a = v/2 the dominating terms are both the C(^*) term in the two-point function, and 
the dominating terms (for a > z//2) in all n-point functions. The distribution becomes: 



e 



-t\{z,)/2 



1 - it/l3c 



□ 



5.2 k G N^\{0, 1}'' density fluctuations 

Now we turn our attention to the density fluctuations of the bare particles. The transforma- 
tion to the quasi-particles ( |2.15| ) has hidden the influence of the external field, therefore it 



plays no role in the distribution functions of the quasi-particles (theorems and |5.5|) . Its 
role will become visible again in the study of the density fluctuations of the bare particles. 
First we treat the fluctuations with high modulation k ^ {0, l}*^ ( ^.3|) . 



Theorem 5.6. The k G N'^yjO, 1}^ density fluctuations (^. Sj ) are Gaussian distributed 
and normal (S = 0) in all phases, the distribution functions are given by: 

• in the normal and critical phases, 

y.(F(iV,))(t)=exp(-t2C(^)/4), 



with ((z) as in ( \5. 14 ), z < e ^^'^ in the normal phase, and z = e ^^'^ in the critical 
phase. 
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• in the condensed phase, 

^iF{Nk)m = exp (-tV8 (2C(^*) +Pocoth(/?i/aV2))) . 

Proof 

Normal Phase: In the normal phase fi^ ^ fi < fi^:,we have the following relation between 
a(0) and 6(0): 



a(0) = 6(0) + LV*^(L) 

where 9{L) = jj^zr^, an unimportant constant converging to a finite value. The fc-mode 
density fiuctuations can be written as the sum over the fc-mode quasi-particle density 
fiuctuations (|5.4| ), Fi = -F(A^{.) and a field fiuctuation Fq 



FLiN,) = L-''/^^J2a\p)a{p + k) + a\p + k)aip) 

p 

= L-"/^^ Yl b\p)Kp + + h\p + k)h{p) + V-'^'^OiL)]- {h\ky^ + h{k)e-'"*') 



2 

p 



= F^ + Fo. 

In the normal phase, the fiuctuations are completely dominated by the first term {Fi = 
FiK)), i.e. 

hm u;i(e"(^^+^«) - e"^^) = 0. (5.28) 

L—>oo 

This can be seen using the following estimates, 

|^^(ei*(i^i+i^o) _ gitFi) I < |^4e^«^iFoe^(*-^)(^^+^°)) | 

Jo 

< tdsujL{e''^^Fy^'Y^\ 
Jo 

using the Cauchy-Schwarz inequality in the last line. The expectation value appearing at 
the rhs in the last line, can be evaluated as follows: 



u;,{e-^^F^e-^^) = j: ^u;,{[F„ F^) 

n=0 

^ E^I^4[^i>^o]n)|, (5.29) 



n=0 

where [F,,F^]o = F^ and [F,,FX+i = [Fi, [F,,F^] 
We now proceed with an estimate of this series. The first term reads 

ul{F^) = L^^-''ujL{h\k)h\k)e-^''^ + h\k)h{k) + h{k)h\k) + h{k)h{k)Q^''^) 
= L^^-'uL{2b\k)b{k) + 1) . 



In I • 
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In order to calculate the higher order terms, remark that we need to keep track only 
of the part h\k)h{k) + b{k)b^{k) = 2b^{k)b{k) + 1, the other terms in Fq and higher 
order commutators of them with Fi contain always an unequal number of creation and 
annihilation operators. Hence, their expectation value is zero and they can be forgotten. 
All contributing terms are of the form b^{q)b{p) and commutation of such a term with Fi 
leads to 

[F^, b^{q)b{p)] = {b\q + k)b{p) + b\q- k)b{p) - b\q)b{p + k) - b\q)b{p - k)) , 



i.e. we find again (at most) 4 terms of the same structure. Using now formula ( |2.18| ) 

\uoL{b\q)b{p)) I < cuL(6t (0)6(0)) , Vp, g G 
we can estimate the rfi^ order term by 

\ul{[F,,FX) I < 2L2^-'^i^L-"'^/24"o.46t (0)5(0)) . 

yielding that the whole series is estimated by 



oo 



^^|u;4[F„Fo\)| < L^^-(2 5^^L— /22"c^46t(o)5(0))+l 

< L^^"" (2a;46t (0)6(0)) e'^^""'" + l) . (5.30) 



n=0 \ n=0 



The term in brackets is bounded for sufficiently large L. Since l?"''^ vanishes to zero, the 
whole term disappears in the thermodynamic limit. This proves ( p.28|) , and the distribution 
of the fluctuations is completely determined by the first term. Hence the distribution of 
F(A^fc) coincides with the distribution of Fi = F{N'^ (|5l^ ), 

lim cu^e'*^^^'^)) = exp(-t2((2)/4)) 
with C(^) as in (|5]T|) and < ;z < e"^'""". 

Critical Phase: A similar argument as in the normal phase can be developed in the 
critical phase {^l = — cL"", with < a < z//2 — 7 and c > 0, some unimportant 
constant) in order to prove that the Fi = F{N'i^) term dominates the fiuctuations. 

The relation between a(0) and 6(0) ( P^.15D depends now on the scaling exponent a (^.9]), it 
reads 

a(0) = 6(0) + -L^+°. 
c 

This yields the following expression for the density fiuctuations 



FL{Nk) = L-'^/2^5^at(p)a(p + A;) + at(p + A;)a(p) 

p 

= L-'/^l J2 h\p)b{p + k) + b\p + k)b{p) + V+''-'''^^{b\ky^ + 6(A;)e-''^) 

p 



Fi + Fo. 
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The Fi term outrules the Fq contribution, i.e. we prove the same formula as ( p.28| ) in the 
critical phase, consider the following estimate 



(t-s)(Fi+Fo) 



^^m+Fo) _ ^tF,^ \ < [ds \ujL{e''^^Foe' 

Jo 

< fdsuL{e''^^F^e-''^^f" 
Jo 



) 



The expectation value ^^(e'^'^^FQ^e ^^^^) can be expanded as 
UL{e''^^F^e-''^^) = ul{F^) +isuL{[F,,F^]) 

- / dsi rds2 ooLie''^^^ [Fu [Fu F^]] e"^^^^^) (5.31) 



The first term on the rhs of (|5.31|) is bounded as 



u;l{F^) = L"+^-^/'ujL{b{kp{k) + b^{k)b{k)) 



£j2a+2'y—i/ ^2 



_/3(,^(fc)-e_j^(0)+cL— ) 



+ 1, 



'l _ e-/3(^L(fc)-eL(0)+cL-°) 

and this goes to zero in the thermodynamic limit since ujL(b^ {k)b{k)^ is bounded if k ^ 
{0, 1}'' and 2a + 27 - z/ < 0, cfr. 

The second term in ( |5.31|) is zero by gauge invariance of the states ujl. 

The third term in ( ^.31j ) can be estimated by a similar procedure as was used in the normal 
phase estimating co'L(e^'^^^FQe~^'^^^), cfr. ( |5.29| )-( |5.3CI| ). The final results is 



|^^(gii(Fi+Fo) _ git(F,)) I < ^ds (o(L2-+27~-) + + 0(L2^+=^"-2'^)e2^^"''') 



1/2 



L^oo 



since 2a + 27 — z/ < and a < u. Hence, ( |5.2(j| ) holds also in the critical phase and the 
distribution of F{Nk) coincides with the distribution of Fi = F^N^.) ( ^.13|) , 

lim u;4e**^(^^)) = eM-t\iz,) /^)), 
with C(-2*) as in (|5.14| ) and z^: = e~^'"^^ . 

Condensed Phase: In the condensed phase the relation between a(0) and 6(0) is: a(0) = 
6(0) + L^/'^^fp^, this yields the following for the fluctuation F{Nk) 

F^Nk) = L--/2^5^at(p)a(p + A;) + at(p + A;)a(p) 

p 

= L-^/^l b\p)KP + k) + b\p + k)b{p) + ^l{b\k)e"^ + b{k)e-'^) 



Fi + Fo 
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In this regime (condensed phase and k ^ {0, 1}''), both terms are equally contributing. 
The fluctuations are normal {6 = 0) and the variances satisfy 

< lim ujl(F?) < oo 
< lim ujlIFq) < oo. 



An argument as in the proof of theorem pT3 can be used here in order to prove that 



lim u;4e'*(^i+^«) = lim cj^e'*^^) lim u}L(e''^°) • (5.32) 

L—>oo L^oo L—>oo 



We use the expansion of the characteristic functions in terms of truncated correlation 
functions ([4.6|) , and prove that all truncated functions vanish 

ujl{F, + Fo, Fi + Fo, . . . , Fi + Fo)^ < 0(L"-^), 

except for the two-point truncated function. 

Odd truncated functions vanish. First note that the one-point truncated function 
vanishes. Consider then the 2n + 1-point truncated function, and suppose that all odd 
m-point truncated functions with m < 2n + 1 are zero. The expression for the 2n + 1 
truncated function reads then: 

lol{Fi + Fo, Fi + Fo, . . . , Fi + Fq)^ = u;l{{Fi + Fo)2"+^) (5.33) 

Use now, 



Fi = ^L~''/y2j2bHp)Kp + k)^ + (^L~''/y2j2b\p + k)b{p)^ 
Fo = {^ob\ky'*>/2) + {^ob{k)e-"f> /2) . 



(5.34) 



The expansion of ^^((Fi -|- Fq)^"'^^) contains terms, leaving out some constants, of the form 

]^-u(r+s)/2 UJL{b\pi)b{pi + k)---b\k)) (5.35) 



pi,... ,pr + s 



where we have r factors of the form b'^{p)b{p + k), s factors 6^(p -|- k)b{p), t factors W{k) 
and u factors &(fc), with r,s,t,u = 1,2,... , 2n -|- 1 and r -|- s -|- t -|- m = 2?t, -|- 1. If we 
calculate the sum of the indices of the creation operators minus the sum of the indices of 
the annihilation operators, i.e. rk — sk + mk — nk, and if r, s, m, n = 0, 1, . . . , 2n -|- 1 and 
r + s + n + m = 2n + l, then 

rk — sk + mk — nk ^ 0. 



Since this is different from zero, it implies that the monome ( |5.35| ) vanishes. If we de- 



compose it into products of two-point functions, all terms contain at least one two-point 
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function with unequal indices for creation and annihilation operator, which is zero, hence 
the whole expression vanishes. By induction, all odd truncated functions ( ^.33| ) vanish. 



All 2n-point truncated functions, with 2r2 > 2 vanish. Consider first the non van- 
ishing term, the two-point function 

u;l{F, + Fo,F, + Fo)j, = ^4(Fi + Fo)2) 

= uj,{Fl) + uj,{F,^) 

= ujl{Fi, + ujl{Fo, Fo)t , 

where the cross-terms ujl{FiFo) ,ujl{FqFi) vanish since the number of creation and annihi- 
lation operators is not equal in those monomes. 

Now consider the 2?T,-point truncated function, 2?t, > 2, and suppose that all 2m-point 
functions, with 2 < 2m < 2n are vanishing. By this assumption we write the 2n-pont 
truncated function as 

uJiiF, + Fo, Fi + Fo, . . .)y = ujl{{F^ + Fo)^") - c,{2n)ujLiF, + Fq, F, + Fq)^ + 0{L"~'') 



where C2{2n) = (2n)!/(2"'n!) is the number of pair-partitions of a set of 2n elements. The 
first term can be expanded as follows. We use first the relation F{Nk) = Fi + Fq and we 
write ujl{F{Ni^Y^) as a sum over terms with the following structure, 

LULiF,F^FoF,Fo---), (5.36) 

where we have 2n factors either being Fq or Fi. Using simple combinatorics, we see that we 
have in total 2^" terms, and (2n)!/(2n — /)!/! terms with < / < 2n factors Fq. We make 
now a further expansion of such an expression using decompositions as in (|5.34|) , we get 
terms consisting of r factors with L~'^/'^ h\p)h{p + k), s factors L~^^'^ h\q + k)b{q), 
t factors b^{k) and u factors b{k), with r,s,t,u = 0, 1, . . . ,2n and r + s = 2n — I and 
t + u = I, i.e. we have terms of the form 

2-2"(po)'/'e'<^(*-")/2L-'^(2"-')/2 J2 Yl ^L{h\k)h\k)h\q^ + k)h{qi)h\pi)h{p^ + k) ■ ■ ■) . 

pi,...,Pr qi,--,qs 



Using similar arguments as in the case of odd truncated functions, we see that the only 
terms which are possibly non-zero are those where r = s = n — 1/2 and t = u = 1/2. These 
terms can be written as a sum of terms consisting of products of two-point functions. The 
sum runs over all pair-partitions of the 4n — / operators in the monome. Since k ^ {0, l}'^, 
such a products of two-point functions contain at most {n — l/2)/2 diverging factors (|5.21| ), 
they are of order 0(L^"/^~'^'"^"~'/^'') and since v > a, such a terms are subdominant. The 
highest order terms are those with the highest number of independent summation indices. 
The non-zero terms with the highest number of independent summation indices have r = 
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{2n — 0/2 summation indices left, they can be constructed if one takes combinations of 
the operators b^{k) and b{k) in factors uji,(b^ {k)b{k)^ or u i,(b{k)b^ (k)^ , and if one combines 
the operators in groups with the same summation index pi of the form b^{pi)b{pi + k) into 
pairs with elements from a group of operators of the form + k)b{pj) and vice versa. 
This can be done for every group of operators since the number of groups of both types 
are equal. It yields factors of the form 

f+{Pi) = ^L{h\pi + k)b{pj + k)) ujL{b{pi)b\pj)) 6{pi - pj), or 
f~iPi) = ^L{b\pi)b{pj)) ujL{b{pi + k)b\pj + k)) 6{pi - pj). 

Terms with less summation indices lead to correction terms which are at most of order 
OIL'"), i.e. the highest order terms are 

(L-72)'- • • • f±iPr)Mb^ik)bik)) ujL[b{k)b\k)) . . . , (5.38) 

where, depending on the order of the factors, we have more or less factors ujL{b{k)b'^ {k)^ 
or u}L{b'^{k)b{ky), and f+{p) or f-{p) ( ^.37|) . But by symmetry in the expansion (|5.34| ), 



there are as many terms in the whole sum starting with f+{pi) as with f-.{pi) or with 
ujL{b^ {k)b{k)) and UL{b{k)b'^ (k)). This holds for all n factors, so in the end we have a 
number of terms which can be written as 

^ (^^(p^) ^ /2 . . . {UiPn-1/2) + f-iPn-1/2)) /2 

pi,--- ,Pn-l/2 

{{uJL{b\k)b{k))+UJL{b{k)b\k))) 12)'" 

and this is equal to 

u,{Fl/2y-"'u,{F^j2f\ (5.39) 
where I is an even number between zero and 2n. 

How many highest order terms ( |5.39| ) are there in the expression for ujl{{Fi + -Fq)^") ? 



There where (2n)!/(2n — /)!/! terms with {] <l <2n factors Fq, cfr. ( |5.36|) , those terms can 
be written as (2n)!/((Z/2)!(n — //2)!)^ terms which where non zero with r = s = n — 1/2 
and t = u = 1/2. These terms are a sum of terms consisting of products of pair-partitions. 
Each original term yields (n — Z/2)!(//2)! terms of leading order ( p.38| ). The total number 
of leading order terms amounts to (2n)!/((//2)!(ra — //2)!) terms of the form (|5.39| ), i.e. 



^L(F(iV,)2") = ^(2n)!/((//2)!(n - I /2)\u;l{FI /2Y-"^ ujL{F'j2f^ + 0(L"--) 



1=0 

(2n)!/2"n! {ujl{F^) + cu^Fq^))" + 0(L°-^) 
C2i2n)ujL{FiNkyy + OiL^-n. 
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Hence, the 2n-point correlation function vanishes in the thermodynamic hmit. The de- 
composition ( ^.32|) is vahd and the hmiting distributions can be written as, 



hm a;4e^*^(^^)) = lim coAe'''^') coAe'''^^) 



2im^e'*^(^^))^e'*^^(^^)) 

exp (-tVS (2C(2*) + Po coth(/3z/(TV2))) . 



□ 



5.3 k G {0, 1}''\{0} density fluctuations 

We continue now with the study of the low-lying k-mode fluctuations. 

Theorem 5.7. The k G {0, 1}^\{0} density fluctuations are Gaussian and normal in the 
normal phase, the distribution function is given by: 

v{F{m)m = exp{-t'a^)/4), z < e-^^-^ 

with ({z) as in ( \5.14 )- 



Proof 

The same argument as in the case of the k G N^\{0, l}'^ density fluctuations (theorem |5 
applies here. 



□ 



Theorem 5.8. In the critical and condensed phases, i.e. when fii = 6^(0) — cL~° with 
c > and < a < z//2 — 7, the distribution of the k G {0, 1}''\{0} density fluctuations are 
different in these regions 

(1) Gaussian and normal (5 = Q) 

^(F(iV,))(t) = exp(-t2c(^,)/4), 



with C(^*) (^s in (\5.14 ), if a < vjl and So; < —27 + v, 
(2) Non- Gaussian and normal (6 = 0) 

72 



. + \2l3c) 



with a{k) = 2(^-'='), if a = z//2 and a + 27 < 0, 
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(3) Non-Gaussian and abnormal (5 = a ~ z^/2j 

/ 1 \ "^^^ 



if a > u/2 and a + 27 < 0, 
(4) Gaussian and abnormal (5 = 7 + 3a;/2 — v/2) 



ip{F{Nk))it) = exp(-t 



if 3a > —27 + u and a + 27 > 0, this last regime includes the the condensed phase, 
i.e. a = = z//2 — 7. 

Proof 

In the critical and condensed phase (/i/, = £^(0) — cL~", < a < a*, c > 0), the relation 
between a(0) and 6(0) reads a(0) = 6(0) + ^L^+". This yields the following expression for 



p 

= L-"/^-^^ J2 h\p)h{p + k) + h\p + k)b{p) + L^+"-/2-5 A(5t(^)e''^ + 6(fc)e-'^) 
p 

= Fi + Fq. (5.40) 



Note that in the condensed phase {a = a^,), we can express Fq as a function of the 
condensate density, i.e. we can substitute ^/po for h/c. 

We have to make a distinction between different regions in parameterspace (|l])-(|D, cfr. 
(Fig. I). 

Region (g) a < iy/2 and 3a < -27 + u 

In this region the Fi term (|5.40|) dominates and the fluctuations can be scaled normally, 
i.e. 6 = 0. The proof is obtained by similar arguments as in theorem ^ for the case of the 
normal and critical phases, we have the estimate 

Jo 

< fdsiJL{e''^'F^e-''^'Y^\ 
Jo 



35 




Figure 4: Representation of different regions in the {a, 7)-parameterspace 
The expectation value on the rhs, can be evaluated as follows 

n=0 

^ E^l^^t^i'^oUl- (5.41) 



w. 

n=0 



The last terms of the series are now estimated as before, using the following estimate for 
the n**^ order term (see 



\ujl{[Fu F^]n) I < /iV2c2L2^-^+2"2"L-"'^/2^46t(o)5(o)) . (5.42) 
The series is bounded by 



°^ I In / °° I 

\S\ I / r „ \ I , ^ , 9 / . O I ^ \ ^ S 



n=0 ■ \ ra=0 

< 0(L2^-'^+3V^"""'). (5.43) 

Since in region (|I]) 3a + 27 — z/ < 0, the expectation value to F^e'^'^^^) is vanishing. 
Hence, the density fluctuations are completely determined by the first term Fi = F{N'i^) 
or: 

hm u;l( e"^^'^^) - e"^(^^)) = 0, (5.44) 
yielding ^{F{N,m = y,(F(iV^))(t), cfr. ( CT) - 
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Regions ©-(H) a > z//2 and a + 27 < 

As in region (|l]), tlie Fi term dominates the Fq term in (|5.4(]| ), the same argument as used 
in region (|l]) can be used here to prove the equahty ( 5.44 ), the only difference being that 
other scahng exponents compared to expression ( ^.42| ) do appear, i.e. 

\ul{[F,,FX) I < /iV2c'^'^2"X-""cu46t (0)6(0)) , 

substituting these exponents in (|5.43|) yields 



5^^M[F„Fo\)|<0(L^--^V^-"). 

n=0 

Since a + 27 < in these regions, this term vanishes in the thermodynamic hmit. Hence, 
also in regions (@)-(0) we have that 

cpiFiN,m = cpiFiN',m, 

and the distribution ip{F{N'f,)) is given by ( |5.15| ) if a = z//2, (region ^ or by ( [5.16| ) if 
a > Z//2 (region ^) . 

Region (g) 3a > -27 + u and a + 27 > 

In this region the Fq term in ( p.40| ) dominates. The variances of Fq are finite, if the scaling 
exponent is chosen as 5 = 7 + 3q;/2 — z//2, while the variance of the first term then vanishes. 
The distribution is completely dominated by the second term, i.e. : 

lim |,;4e^*(^^+^")-e'*(^"))|=0. 

L^oo 

This is proved using similar bounds as before: 

Jo 

< fdsu;L{e''^'>Fy^'^y^\ 
Jo 

The expectation value appearing is again expanded as: 

u;,{e^^^^F^e-^^) = J2 ^u;,{[Fo, F^n) • 

n=0 

An easy calculation learns that this series is cut off after the third term, 

^4e'^^»F2e-'^^«) = u^iF^) + u;l{[Fo, F^]) + cu^i^o, [Fq, F^]]) (5.45) 
The first term reads 

u;l{F^) = L"'-^'J2uL{b\p)b{p))+2u;L{b^{p)b{p))uL{b\p + k)b{p + k)) 



V 



+ UL{b\p + k)b{p + k)) 

pe{o,iY p^lo,!}"" 
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The terms in the sum over p G {0, lY are at most of order 

UL{bHp)b{p)) UL{h\p + mp + k)) = 0(1^'^), 

and since 27 + 3a > 2a in region they are vanishing; the second term is a Riemann 
sum converging to a finite integral, 

hmL-'^ ^ ... = C(^*), 

with C(-2*) as in ( p.l4| ). The scahng exponent satisfies z/ — 27 — 3a > in region (^), 
therefore this term vanishes in the thermodynamic hmit as welL 

The second term in (|5.45|) is zero because the constituents are expectation values of 
monomes with an unequal number of creation and annihilation operators 

iuL{[Fo,F^]) = iujL{F^[Fo,F,] + [Fo,F^]F,) 

= AL-^-2"cuL(iFi (6(2fc)e-^<^ - b^{2k)e''^ 
+ 6(0)e-''^-6t(0)e''^) + h.c.) 

= 0. 



The third term of ( 5.45|) reads 



ul{[Fo, [Fo, F^]]) = ujl{2[Fo, F,]^ + F,[Fo, [Fq, F,]] + [Fq, [Fq, F,]]F,) . 

The expectation values of the second and third term in this expression vanish since 
[Fo, [Fo, Fi]] = 0, the first term is expanded as 

u;l{[Fo,F,]') = -L-^"-2^^(a;46t(o)6(o)) + l+c.46t(2A;)6(2A;))) 
and tend to zero in region (^. 

The higher order terms in ( ^.451 ) vanish since they all contain the factors [Fq, [Fq, Fi]] = 0. 
Hence the distribution of these density fluctuations coincides with the distribution of the 
field fluctuations, i.e. Fq = ^F(A+) cfr. ( CT) , 



^{F{N,)m = hm ^4e^*^(^^)) = lim u.fe^^^'^^^A = exp (-t'-^] 



□ 
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5.4 Unmodulated {k = 0) density fluctuations 
Theorem 5.9. The limiting distributions of the k = density fluctuations ( [5. j| j are 
• Gaussian and normal 

^{F{No)m=expi-t'az)/2) 



with ({z) as in ( \5.14 ), in the normal phase and in the critical or condensed phases 



(with /iL = eiiO) — cL ") if a < iy/2 and if 27 + 3a; < u. 
Non- Gaussian and normal 



V.(F(iVo))(t)=e-*'^(-)/2 



1 - it/f3c 



with C(-2*) fls in ( |5. 1^ ), in the critical phase if a = 1^/2 and 7 < — z//4, 



Non- Gaussian and abnormal (5 = a — v /2) 

^{F{No)m = 



1 - it/l3c 

in the critical phase if a > vj2 and a < —27. 
Gaussian and abnormal (6 = j + 3a/2 — i'/2) 



^(F(iVo))(t)=exp(-t^^) 



in the critical or condensed phase if a > —27 and a > u/S — 2'-/ /3. In the condensed 
phase, i.e. if a = a* = z//2— 7, this function can be written in terms of the condensate 
density po, i.e. substitute p^J'^ / I3h for h? j fic' . 



Proof 



The proof goes along the same hnes as the proofs of theorem 5. / and theorem 5.8 



□ 



6 Concluding Remarks 

Clearly our key contribution consists in the rigorous analysis of the field and density fiuctu- 
ations in the Bose gas with attractive boundary conditions. Explicit distribution functions 
of the fiuctuations are computed which are as well of the Gaussian or the non-Gaussian 
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type, of the normal as well of the abnormal type. These results, clearly indicate the influ- 
ence on the details of the boundary conditions, like the external field ( p.l3|) , the elasticity 
of the boundaries ( p.2| )-( p3D , and the volume dependence of the chemical potentials (|3^). 

In the technique to study the thermodynamic limit we adopt here, the interplay between 
the scaling of the external field ( p.l3| ) and the convergence rate of the chemical potentials 
( p.8|) determines the limiting thermodynamic phase. This interplay also determines the 
distribution and the scaling of the field and density fluctuations. For the field fluctuations 
(theorem ^]2| ) and for high- modulated density fluctuations (theorem |5l6|) , the structure 
is determined by the division between normal, critical and condensed phases. For low- 
modulated and unmodulated density fluctuations (theorems ^]8|and p.9| ), different regions 



appear in the critical phase where these density fluctuations have different scaling laws and 
distribution functions (cfr. Fig. ^). Hence, the study of the Bose gas with external field, 
cannot be confined to Dirichlet boundary conditions or a 7 = scaling of the external 
field. As our explicit calculations demonstrate, going beyond these limitations reveals 
an unexpected richness in the fluctuation distributions. This suggests that a choice of a 
particular strength of the external field should be well motivated on physical grounds. 

The results about the (unmodulated) density fluctuations can be compared with previously 
obtained results P] . There is a correspondence on a heuristic level in the sense that density 
fluctuations are Gaussian and normal in the normal phase, non-Gaussian and abnormal in 
(part of) the critical phase and Gaussian but abnormal in the condensed phase, and the 
form of the explicitly obtained non-Gaussian distribution ( p.26[ ) is of the same structure as 
the ones calculated by Angelescu et. al. in [Q. But the scaling exponents they calculated 
are different. They found, in three dimensions respectively 6 = 0, 0.5, 3.5 for the normal, 
critical and condensed phases, whereas we found respectively 5 = 0, 5 G [0,1.5), and 
S G (0,1.5), indicating that a difference in the boundary conditions changes the scaling 
behaviour as well. Furthermore, we found that density fluctuations can also be Gaussian 
or normal or both in the critical regime, depending on different choices of parameters. 
Sufficiently strongly modulated fluctuations, i.e. if k ^ {0, lY are always Gaussian and 
normal, confirming that there is a substantial Gaussian element even at or below the critical 



point, cfr. the discussion about the classical Curie- Weiss model in |T3|,|T^]. 

The explicit form of the distribution functions can be obtained because of the quasi-free 
character of the equilibrium states. This model is original in the sense that the spectrum 
shows an energy gap in the case of attractive boundary conditions. Previous explicit studies 
of fluctuations where only performed in the case of spectra without an energy-gap |ll|-10 



One would also like to know whether such changes in boundary conditions are as important 
in models for interacting Bose gases. 

Finally, we want to remark that there are more open problems to be studied in this model. 
In particular we considered here a special thermodynamic limit using a power low depen- 
dence of the chemical potential on the volume, yielding explicit dependence of the critical 
exponents on this power law. An other natural way of taking the thermodynamic limit is a 
limit by taking the density constant. It is a question whether in this case the different dis- 
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tribution functions can be computed explicitly and whether their behaviour is analogous. 
Our experience by now is that properties of fluctuations are very dependent on the type of 
thermodynamic limit taken. This problem is related to the question about the equivalence 
of ensembles for the free gas with a finite gap in the energy spectrum, this is currently 
under investigation, and hopefully sheds new light on this problem. 
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